Genetic tracing of market wildlife and viruses at the epicenter of the COVID-19 pandemic (preprint)

Zoonotic spillovers of viruses have occurred through the animal trade worldwide. The start of the COVID-19 pandemic was traced epidemiologically to the Huanan Wholesale Seafood Market, the site with the most reported wildlife vendors in the city of Wuhan, China. Here, we analyze publicly available qPCR and sequencing data from environmental samples collected in the Huanan market in early 2020. We demonstrate that the SARS-CoV-2 genetic diversity linked to this market is consistent with market emergence, and find increased SARS-CoV-2 positivity near and within a particular wildlife stall. We identify wildlife DNA in all SARS-CoV-2 positive samples from this stall. This includes species such as civets, bamboo rats, porcupines, hedgehogs, and one species, raccoon dogs, known to be capable of SARS-CoV-2 transmission. We also detect other animal viruses that infect raccoon dogs, civets, and bamboo rats. Combining metagenomic and phylogenetic approaches, we recover genotypes of market animals and compare them to those from other markets. This analysis provides the genetic basis for a short list of potential intermediate hosts of SARS-CoV-2 to prioritize for retrospective serological testing and viral sampling.

Tracking SARS-COV-2 Variants Using Nanopore Sequencing in Ukraine in Summer 2021 (preprint)

Since spring 2020, Ukraine has experienced at least two COVID-19 waves and has just entered a third wave in autumn 2021. The use of real-time genomic epidemiology has enabled the tracking of SARS-CoV-2 circulation patterns worldwide, thus informing evidence-based public health decision making, including implementation of travel restrictions and vaccine rollout strategies. However, insufficient capacity for local genetic sequencing in Ukraine and other Lower and Middle-Income countries limit opportunities for similar analyses. Herein, we report local sequencing of 24 SARS-CoV-2 genomes from patient samples collected in Kyiv in July 2021 using Oxford Nanopore MinION technology. Together with other published Ukrainian SARS-COV-2 genomes sequenced mostly abroad, our data suggest that the second wave of the epidemic in Ukraine (February-April 2021) was dominated by the Alpha variant of concern (VOC), while the beginning of the third wave has been dominated by the Delta VOC. Furthermore, our phylogeographic analysis revealed that the Delta variant was introduced into Ukraine in summer 2021 from multiple locations worldwide, with most introductions coming from Central and Eastern European countries. This study highlights the need to urgently integrate affordable and easily-scaled pathogen sequencing technologies in locations with less developed genomic infrastructure, in order to support local public health decision making.

Accuracy in near-perfect virus phylogenies (preprint)

A bstract Phylogenetic trees from real-world data often include short edges with very few substitutions per site, which can lead to partially resolved trees and poor accuracy. Theory indicates that the number of sites needed to accurately reconstruct a fully resolved tree grows at a rate proportional to the inverse square of the length of the shortest edge. However, when inferred trees are partially resolved due to short edges, “accuracy” should be defined as the rate of discovering false splits (clades on a rooted tree) relative to the actual number found. Thus, accuracy can be high even if short edges are common. Specifically, in a “near-perfect” parameter space in which trees are large, the tree length ξ (the sum of all edge lengths), is small, and rate variation is minimal, the expected false positive rate is less than ξ /3; the exact value depends on tree shape and sequence length. This expected false positive rate is far below the false negative rate for small ξ and often well below 5% even when some assumptions are relaxed. We show this result analytically for maximum parsimony and explore its extension to maximum likelihood using theory and simulations. For hypothesis testing, we show that measures of split “support” that rely on bootstrap resampling consistently imply weaker support than that implied by the false positive rates in near-perfect trees. The near-perfect parameter space closely fits several empirical studies of human virus diversification during outbreaks and epidemics, including Ebolavirus, Zika virus, and SARS-CoV-2, reflecting low substitution rates relative to high transmission/sampling rates in these viruses.